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ABSTRACT 

We propose a simple model for the origin and evolution of Lya forest clouds based on 
cosmological conductive/cooling fronts. In this model the Lya forest arises in the interfaces 
between the IGM and cold clouds that could be tentatively identified with protogalaxies. 
Most of the properties of the Lya forest absorbers are reproduced with a very restricted 
number of assumptions. Among these are the correct range of HI column density, cloud 
sizes and redshift and HI column density distributions for the absorbers. Several predic- 
tions and implications of the model are briefly discussed. 

Subject Headings: intergalactic medium — quasars: absorption lines 



1. Introduction 



The remarkable observational progresses made in the last few years allow to tackle the 
problem of the origin and evolution of Lja forest clouds in a more firm and quantitative 
manner. In brief, the available observational constraints that every model for the Lya for- 
est clouds should fulfill are the following. The typical range of hydrogen column density 
is 10 12 cm~ 2 < N H i < 10 16 cm -2 (Dobrzycki & Bechtold 1996). Doppler parameters are 
typically b = 30 km s _1 , even if some claims of much smaller values of b have been made 
(Pettini et al. 1990; Carswell et al. 1991; Rauch et al. 1993; Giallongo et al. 1993). Lower 
limits on the sizes are provided by a small number of experiments using quasar pairs (Din- 
shaw et al. 1994, Bechtold et al. 1994) which generally agree that transverse sizes should 
be > 100 kpc. 

Here we propose an alternative view of the origin and evolution of Lya forest absorbers, 
based on the theory of conductive/cooling fronts described by Ferrara & Shchekinov (1993) 
(FS93). The basic idea is that Lyct forest clouds may represent the interface region between 
the IGM and cold clouds that could be tentatively identified with protogalaxies. Since the 
IGM is in a thermally unstable state, the joint action of thermal conduction and radiative 
instabilities drives a cooling wave into the hot IGM. Thus, contrary to the common view 
in which thermal conduction should eventually lead to the evaporation and destruction 
of the cold phase, we argue that protogalaxies initiate the condensation process of gas 
which can be subsequently accreted by the galaxy itself. The relatively rarefied transition 
layer between the two phases represents a suitable environment for the formation of the 
Lya forest forest absorption lines. Although very simple, the model is able to reproduce 



most of the properties of Lya forest clouds with a very restricted number of assumptions. 



2. Conductive/Cooling Fronts 



In this Section we describe briefly some results concerning conductive/cooling (CC) 
fronts following FS93 and Ferrara & Shchekinov (1996). We define a CC front as the 
interface between a hot and cold gas phase, whose structure is governed by the combined 
effects of thermal conduction and radiative cooling. As shown by FS93, dynamical effects 
(i.e., shocks) can be important at the initial stages of the front evolution when the conduc- 
tive timescale is shorter than the dynamical timescale of the system. However, after this 
transient the evolution relaxes to a quasi-isobaric regime. We will therefore consider only 
evolved stages for which pressure equilibrium holds approximately. For the same reason, 
we consider only "classical" Spitzer conductivity since saturation effects are important 
only in the initial evolutionary phases (FS93). We neglect magnetic fields, mostly because 
their existence and strength at high redshift are far from being established. 

In the isobaric approximation, the hydrodynamic equations describing a CC front reduce 
to the energy equation, which could be cast in Lagrangian form as 



where T is the gas temperature, C is the net cooling rate per unit mass, c p is the specific 



coefficient, is the reduced thermal conduction coefficient. The Lagrangian mass variable is 



&T c~ x L 9 \ 9T ] - 
dt Cp dq * dq 



(2.1) 



heat at constant pressure; x = k p/ c p, where k 



rjT 5 / 2 is the classical thermal conduction 



q = / dxp(x, t) 
Jxo 

where for all t the coordinate of a "reference" particle xo(t) is determined by the condition 
v(xo,t) = 0. In eq. (2.1) lengths are normalized to the Field length (Field 1965) (La- 
grangian formulation), q F = (xTc p / C) 1 / 2 , and time to the cooling time t c = c p T '/ both 
calculated for the hot medium. 

The cooling function depends on the details of the microscopic processes responsible for 
the heating and cooling of the gas; for our purposes, however, it will be sufficient to consider 
a general form of £ reproducing a two phase medium with a cold stable, and a hot unstable 
equilibrium. This assumption closely resembles the IGM, supposedly constituted by cold 
clouds in pressure equilibrium with a surrounding, thermally unstable, hot diffuse gas. 
The simplest analytical function that retains such characteristics is C(T) = T(l — T); the 
point T = 1 is thermally unstable, (dC/dT = — 1), while T = represents the thermally 
stable phase. In addition to the equilibrium points, this function also correctly mimicks 
the negative slope of the actual cooling function in the range 10 4 K 10 6 K. 

Fig. 1 shows the temporal evolution of the CC front, as obtained from the numerical 
solution of eq. (2.1) with the above cooling function. Initially a cold cloud at T = is 
immersed in the hot gas at T = 1. The evolution of the CC front induces a thermal wave 
in the cold medium and a cooling wave into the hot medium; the temperature profile is 
very steep due to the strong nonlinearity of the thermal conduction (shown by the enlarged 
view in Fig. 1). At about t ~ t c the thermal wave propagation is inhibited by radiative 
losses (see FS93), and the cloud now acts as a finite perturbation for the unstable hot gas, 



whose only possible fate is to condensate. As the CC front moves away from the origin, a 
cooling wave develops (t ~ 2t c ); the profile of the wave becomes stationary and its shape 
corresponds to a quasi-linear thermal conductivity, in that (dT/dq) 2 <C \T(d 2 T/dq 2 )\. 
This fact suggests to compare the numerical results with an analytical travelling wave 
solution of eq. (2.1) obtained for linear conductivity (x = const.). Making the position 
T{q,i) = T(£ = q + 5Ut), substitution into eq. (2.1) yields the solution 

where U is the propagation speed of the wave. In Fig. 1 we have fitted the curve of 
eq. (2.2) to the cooling wave, adjusting the parameter U to take into account that the 
exact solution is only quasi-linear; the agreement is remarkable for U = 4/3 = e -1 . Of 
course, it is not surprising that the value of U is close to unity, since it can be seen from 
dimensional analysis that it must be U ~ q F jt c = 1, in our units. The evolution depends 
rather weakly on the detailed shape of C as long as an unstable point does exist. We will 
use this approximation to study the interface bewteen Lyct forest clouds and the IGM in 
the next Section. 

3. Conductive/cooling model for Lya forest absorbers 

To model Lya forest clouds as cosmological CC fronts it is necessary to adopt a model 
for the IGM. Ikeuchi & Ostriker (1986) (and more recently Giroux & Shapiro 1995) have 
studied in detail the thermal history of the IGM from the era of the reheating to the present; 
they conclude that the most plausible scenario requires that both photoionization and bulk 



mechanical heating (i.e., shocks) contribute to the heating. The IGM reionization epoch is 
still uncertain, but the limits on the Gunn-Peterson effect (Webb et al. 1992) suggest that 
it must have occurred at z > 5. Ikeuchi & Ostriker show that, if the reionization occurred 
at z = 10 - for example - the IGM entered the adiabatic expansion phase already at z ~ 4. 
Since we will mostly concentrate on the z ^ 4 epoch, we will assume in what follows that 
the IGM is adiabatic, and therefore temperature and pressure are T(z) = T (l + z) 2 ■, 
p(z) = po(l + z) b ; the local values are T = 3 x 10 4 K and po = 3 x 10~ 2 cm -3 K; we use 
the cosmology = 1 and h wo = 0.5; in addition & = 0.03. 

Since the heat transport relies on the presence of a hot IGM, it is likely that CC fronts 
were generated immediately after the reionization epoch. Next, they propagate away from 
the cloud with velocity decreasing rapidly with redshift: U ~ q F /t c oc (1 + z) 6 , assuming 
that IGM cooling is dominated at early epochs by inverse Compton scattering on the 
microwave background. A stationary cooling wave profile, similar to the one given in eq. 
(2.2), will form after a few cooling times, as demonstrated in Fig. 1; this requires that the 
Hubble time is larger than t c , a condition satisfied for z ^ 4. If reionization occurred at 
z = 10, the distance travelled by the cooling wave at z = 4 is d = (m^o)" 1 dz(l + 
z)~ 5 / 2 UrT 1 = 60 kpc; the cool material behind the wave is accreted by the parent 
protogalaxy at a rate Mi n ~ mHnd 3 tJ I 1 4 = 5 x 10 7 M Gyr -1 , where tn,4 is the Hubble 
time at z = 4. As the front has reached the steady-state, its width (of the order of q F ) 
is regulated only by the changes of the IGM parameters due to Hubble expansion; since 
q F oc (1 + z) 2 , the front shrinks as z decreases (it can be shown that the interface reacts 
almost instantaneously to external changes). 



The neutral hydrogen column density of the Lyct forest cloud, N HI , will be in general a 
function both of the impact parameter, 6, at which the interface is intersected by the line 
of sight and of redshift. Assuming spherical symmetry, with radial coordinate r{q), the 
expression for Nhi is 

N HI (b,z) = — [ dq{l-x[T(q,z),n(q,z)]}^^L==. (3.1) 
n>H J \Jr 2 (q) - b z 

The hydrogen ionization fraction x has been derived under the assumption of ionization 
equilibrium due to photo- and collisional ionizations and that the time variation of the 
diffuse UV flux is J(z) = 10~ 21 [(1 + z)/3.5] 2 ergs cm -2 s -1 Hz -1 sr -1 . The temperature 
profile follows from eq. (2.2) and n(q) = p/kT(q); we have neglected geometrical correc- 
tions in the solution T(q) due to spherical geometry. The lower integration limit is the 
Lagrangian impact parameter implicitly defined by mjj 1 dq/n(q) = b — d, while the 
upper limit is the characteristic size of the interface. The size of the central cold core is 
determined by the distance d ~ 60 kpc travelled by the front in the early stages of evo- 
lution; however, the value of the integral depends very weakly on any reasonable choice 
for d. A sketch of the geometry is given in Fig. 2. Impact parameters b < d are not 
considered because of the negligible contribution of the Lyct forest cloud to the central 
damped Lyct system. 

Fig. 3 shows the curves for N HI from eq. (3.1) as a function of q(b) for different 
redshifts. Therefore, q(b) = corresponds in physical units to an impact parameter b = d 
and q(b) = 1 corresponds to b = e£ F + d, where t F = q F /mHn is the Field length in 
physical units; £ F depends on redshift as 



t r = {^) ' = 333 C 1 + z )~ l k P c > ( 3 - 2 ) 

for the adopted IGM parameters. 

The hydrogen column density is a rather flat function of q(b) for values q(b) < 0.1, 
i.e., close to the central cold cloud, and decreases roughly as q(b)~ 4 for z = 4 and less 
steeply for smaller values of z. The values of N HI are in the range 10 12 — 3 x 10 16 cm -2 
for < z < 4. An useful analytical approximation to the curves shown in Fig. 3 between 
z = and z = 3 is 



N H i(q) = N° HI 
with N^j = 7 x 10 12 cm" 2 , and a = 0.1. 



(1 + *) 



(3.3) 



4. Implications 

In the present model Nhi is a function both of redshift and of the impact parameter 
for the cloud: it decreases with decreasing z and with increasing b; the variation range is 
10 12 cm" 2 < N HI <3x 10 16 cm" 2 for < z < 4. Therefore, below Nhi = 10 12 we expect 
a turnoff in the distribution of the clouds; this prediction awaits an observational test from 
oncoming high sensitivity observations with the Keck telescope. 

The temperature is a function of the depth inside the cloud along the line of sight (see 
Fig. 2), with radial dependence as described by eq.(2.2). Its range is bracketed by the 
temperature of the cold cloud (here assumed to be T c = 10 4 K) and the temperature of 
the IGM, T(z) (roughly 10 5 — 7 x 10 5 K). For a given Lyct forest cloud, the temperature 



is anticorrelated with Nhi, since the inner regions are colder than the external ones. 
However, since this cold material is far from hydrostatic equilibrium in the gravitational 
field of the parent protogalaxy, it will be accelerated towards the center acquiring bulk 
velocity that could reverse the sign of the T — Nhi correlation. This conclusion must be 
substantiated by fully hydro dynamical calculations. 

An interesting feature of the model is the presence of a natural scale for the size of the 
clouds, the Field length £ F (eq. [3.2]). Thus, the size of Lyct forest clouds is constrained 
directly by physical arguments. For the adopted model of the IGM, typical transverse sizes 
are ~ 50 kpc and ~ 150 kpc, for z = 4 and z = 1, respectively. Hence, for a constant 
number of absorbers per comoving volume, the probability to detect a Lya forest cloud 
in a quasar pair should be higher at low redshifts if the signal-to-noise ratio were infinite 
(see below). 

The expected redshift distribution of absorbers can be calculated as follows. The number 
of absorbers between z and z + dz is 



where N a = AT (1 + z) 3 is the density of absorbers. The maximum impact parameter, 
b m , is determined by the observational threshold for detection of Lya forest absorption: 
b m = b(N H l I ). For example, from Fig. 3, at z = 3 and for Njjj = 10 13 cm -2 , which 
represents a typical threshold for current observations, the maximum value of q(b) sampled 
is ~ 0.4 (there is an almost linear relation between q(b) and b). For the previous value of 
Njjj, a fit to the numerical results gives b m oc £ F (1 + z) a . In 1 < z < 4 the best fit value 



b 





o 



is a = 1.4. For a higher value of Njjj, the corresponding value of a increases, i.e., the 
distribution becomes steeper. Upon substitution in eq. (4.1), with a = 1.4, we obtain 

dN/dz oc (1 + z) 1 - 3 . 

The index of the distribution is slightly smaller than the one found by Bechtold (1994) 
(7 = 1.32 ± 0.24) for a sample of lines complete to equivalent width Wth = 0.16 A. 
According to her curves of growth, depending on the adopted Doppler parameter, this 
corresponds to 6 x 10 13 cm~ 2 <, Nhi ~ 10 14 cm~ 2 , about an order of magnitude larger than 
Njjj. Since Bechtold has also noted a decreasing trend for 7 with decreasing W t h, it is 
reasonable to expect that for lower W t h the distribution flattens to a value similar to 
the one we have found. However, this physical effect can be difficult to disentangle from 
spurious line blending effects (Trevese et al. 1992). 

The column density distribution for the absorbers is 

3$^ = ™. / dz b(N H ,,z)^(N H ,,z)-jp—(N H i,z)(l+zf^. ,4.2) 



Using the definition of the Lagrangian coordinate q to eliminate b and eq. (3.3) to eliminate 

q, we get dN/dN HI oc iV^J/(JV HJ , z m ). 

Numerical integration shows that the function f(Nui, z m ) depends very weakly on Nhi 
and can be approximated by a powerdaw with index —0.1. Therefore, from eq. (4.3) we 
find that the column density distribution for the Lya forest clouds is also a power-law with 
index j3 = —2.1. The value of j3 is independent of z m and somewhat higher than the one 
(13 = -1.4) obtained by Dobrzycki & Bechtold (1996). 

Two additional consequences of our model can be outlined: i) Lya forest clouds should 



be quite often associated with damped Lya clouds, from which they are originated, at 
least at high z. At low z, damped Lya systems can be destroyed by violent episodes of 
star formation (Salpeter 1995) and isolated "fossile" Lya forest clouds remain; ii) metals 
abundances at the level detected by Tytler (1995) could be easily reproduced by our model 
if part of the Lya forest gas is coming from a metal-rich protogalaxy, or if the IGM is 
polluted by metals associated with the same blast waves responsible for its bulk heating. 
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Figure 1 Evolution of the structure of a conductive/cooling front; time is in units of the 
cooling time, t c , in the hot gas; temperature is normalized to the temperature in the hot 
gas. Solid curves refer to (from right to left) t = 0, 0.5, 1, 2, 3, 10; the dotted curve shows 
the analytical approximation for t = 10. The inner panel shows an enlargement of the 
region close to the origin. 



Figure 2 Sketch of the geometry used for the conductive/cooling model. 



Figure 3 Distribution of the hydrogen column density for a Lya forest cloud as a function 
of the Lagrangian impact parameter q(b) (see text) for different redshifts z = 4,3,2,1,0 
from the uppermost to lowest curves; q(b) is in units of the Lagrangian Field length at the 
relevant redshift. 



